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4^ Abstract 
Oh 

Neutrinoless double beta decay can in general produce electrons of either chirality, in 

r-| contrast with the minimal Standard Model (SM) extension with only the addition of the Weinberg 
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operator, which predicts two left-handed electrons in the final state. We classify the lepton number 
violating (LNV) effective operators with two leptons of either chirality but no quarks, ordered 
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according to the magnitude of their contribution to Oi//3/3 decay. We point out that, for each 
5 o, .e ..ee .... a., ... ... . o„. ..V ope... . 

^ corresponding type to lowest order, and these have dimensions 5, 7 and 9, respectively. Neutrino 



masses are always induced by these extra operators but can be delayed to one or two loops, 
depending on the number of RH leptons entering in the operator. Then, the comparison of the 
decay rate and neutrino masses should indicate the effective scenario at work, which confronted 
with the LHC searches should also eventually decide on the specific model elected by nature. We 
also list the SM additions generating these operators upon integration of the heavy modes, and 
discuss simple realistic examples of renormalizable theories for each case. 
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I. INTRODUCTION 



The remarkable observation of neutrino oscillations (see [1] and [2, 3] for recent reviews) 
provided the first direct evidence of physics beyond the Standard Model (SM), these effects 
are best explained by providing the neutrinos by small masses and appropriate mixing angles 
[4, 5] (see also [6] for a recent fit). In contrast with the quark sector, however, neutrino 
masses are not necessarily of the Dirac type, yet oscillation experiments are not sensitive 
to determine whether neutrino masses are of the Majorana type [7]. Fortunately there 
is another process, neutrino-less double-beta (Oz//3/3) decay ([8, 9] and [10] for a review) 
that probes this property of the neutrino sector, and has achieved sufficient sensitivity 
(see [11-13] for recent reviews) to provide interesting constraints on the lepton-number 
violating (LNV) processes that can produce it. In this paper we will be concerned with 
general properties of both neutrino masses and Oi//3/3 decay; we will strive to provide model- 
independent description of the effects that concern us, yet we will make connection with 
specific models that provide concrete and important illustrations of the arguments presented. 

Since both neutrino masses and Oz//3/3 decay are low energy processes, an effective La- 
grangian approach [14-19] is the proper starting point of any model-independent discussion; 
in it, all virtual new physics (NP) effects are parameterized by the coefficients, C*^"-', of the 
corresponding effective operators O*^"^ Explicitly, 



where n denotes the canonical dimension of the operator, i labels the independent operators 
of a given order and A the NP scale. The O^-"'* respect all the local symmetries of the SM, 
but not necessarily the global ones; when the NP contains several scales A > A' > ■ ■ ■ the 
coefficients C may contain powers of A'/A (we will assume that the new physics is weakly 
coupled and decoupling). 

Several earlier papers [20-23] have followed this approach, considering, however, only 
effective interactions that do not involve the SM gauge bosons. Here we consider a different 
class of theories where the NP does not couple directly to the quark sector ^, so that the 
effective interactions involve only leptons and gauge bosons (coupling to the gauge bosons is 

^ Operators involving leptons and quarks with no gauge bosons generate neutrino masses at 1 to 4 loops [23- 
26] and may receive enhancements from top Yukawa couplings. 
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generated whenever the NP is not a SM gauge singlet). We will show that in this case one can 
provide a simple classification of the effects that concern us in terms of only three operators, 
each of which can be generated at tree level by different types of NP. In the unitary gauge 
these operators give the vertices i^li^l, WcrVl and WWcrCr and have dimension 5,7 and 
9, respectively. 

This allows for three scenarios wherein one of the operators is generated at tree level 
and the others via loops. A simple example exhaustively considered in the literature has 
tree-level generated neutrino masses via a high-scale see-saw mechanism [27-30], with ef- 
fective WcrUl and WWeRCR vertices generated radiatively. Here we concentrate on the 
complementary situations where the other operators are generated at tree level and neu- 
trino masses are generated radiatively. Though we will discuss the three cases separately it 
is of course possible for the NP to simultaneously contain all of them. Indeed, in specific 
models, phenomenological constraints might necessitate such complications; this is the case 
in realistic left-right (L-R) models [31]. However, the interesting result of the analysis we 
present is that, whatever the complications, they reduce to combinations of the three cases 
discussed below whenever our assumptions are applicable. 

The classification of NP contributions to Oz//3/3 decay is straightforward. The final state 
involves two electrons of either chirality which may proceed directly from the effective vertex 
generated by the NP, or from the WvlCl SM vertex. As a result LNV effects contributing to 
Oi//3/3 decay can be classified according to the chirality of the final leptons. It is remarkable 
that, as we shall show, there is only one low est- order operator for each of the three possible 
chirality assignments exhibiting a notable connection between chirality and dimension: ^ 

(9(5) = for two LH leptons (LL), (2) 

C)(^) = (0tD'^0)(0te^^^4), for one LH lepton and one RH charged lepton (LR), (3) 
(9(9) = e]ie%{(j)^D^'4)){(f)^Df,4)), for two RH charged leptons (RR). (4) 

Here we omitted fiavor indices and denote the light scalar isodoublets by 0, the left-handed 
(LH) lepton isodoublets by i^, and right-handed (RH) lepton isosinglets by e^; we also use 
= i(J2<P* 1 and ^L = i(J2^'L- The electric charge is the sum of the third component of isospin 

^ Notice that the O^^^ lepton number (LN) assignment is opposite to that of O'^'^^ and O*^^^ in order to 
preserve the usual convention for 0^^\ 
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FIG. 1: One (left) and two (right) loop neutrino masses for LR and RR operators, respectively. 
Arrows indicate fermion number flow. 

and the hypercharge, Q = 1^ + Y ; thus Cl and have hypercharge —1/2, while and 
have hypercharge +1/2. 

The operators in Eqs. (2-4) also provide contributions to the neutrino masses; the first 
one at tree level, as first noted by Weinberg [32] (see also [33]), the other two, radiatively at 
one and two loops, respectively; see Fig. 1. We will show that ^ (here a, h are family indices) 

{m^)ab oc ^C^?*, for LL, (5) 

{'^^)ab oc {ma&2 + "^bC^Ha) > LR, (6) 

{m^)ab oc ^^gJ^^^^ m^C^fc^mfe, for RR, (7) 

where the vacuum expectation value (VEV) (0) = f ~ 174 GeV is the electroweak symmetry 
breaking parameter and b the corresponding charged lepton masses. It is worth noting 
that the contributions from (9*^^'^^ to niy have a natural hierarchy derived from the charged 
lepton mass factors and may be used to generate textures naturally. We will discuss this 
issue in a forthcoming publication. 

In all cases the neutrino masses vanish when the scale of NP A — )■ oo, as required by the 
decoupling theorem [34]. This implies that the neutrino oscillation data imposes a upper 
bound on A, while the limit from Oz//3/3 decay provides a lower bound on this scale. It 
also follows that the NP cannot be flavor blind, so that specific models in general will be 

^ The dominant contributions obtained in specific models are proportional to A^^ since the lowest dimension 
operator describing neutrino masses is 0'^^\ and this operator has dimension 5. Once LN is broken, O^^^ 
will always be generated, at one or two loops (barring model dependent cancellations), even if it does not 
appear at tree level; this happens in the models we will consider that generate O^^-' or O^^' at tree level. 
Deriving these estimates requires care; in particular, estimates obtained using the unitary gauge are often 
not reliable. 
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constrained, not only by the neutrino and Ou/S/S data, but also by lepton-flavor violation 
(LFV) constraints [35]. 

If the NP generates any one of these operators at tree level, the remaining two will be 
generated at one or two loops; when O^^^ is generated radiatively the neutrino masses are 
proportional to the coefficients of the tree-level generated operator and are, in this sense, 
predictable. In a companion paper [35] we provide a realistic and highly constrained model 
that illustrates this scenario, having only O^^^ generated at tree level while 0^'^\ O^^^ appear, 
respectively at one and two loops. In particular, neutrino masses are generated at two loops 
(see Eq. (7)); the wide literature on radiative neutrino masses provides many additional 
examples (see for instance [36-39]). 

We will see that current data on Ou/SP decay implies that the scales associ- 
ated with NP generating 0^^\ O^''^ and O^^^ at tree level are, respectively A > 
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C, 



(5) 



TeV, 102 



a 



(7) 



1/3 

TeV, and 
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(9) 



1/5 

TeV. For first type of NP we expect all 



collider effects to be negligible; while we expect that the NP responsible for O^^^ will be 
probed at the LHC [35]. The intermediate case of O(^) may or may not have collider signa- 
tures, depending on the details of the model; we examine one such case below. 

In next section we shall define our notation and using the effective Lagrangian approach 
we classify the lowest order interactions (operators) invariant under the SM gauge group 
with two external leptons but no quarks, mediating Oz//3/3 decay (it turns out that such 
operators have equal number of covariant derivatives and RH leptons). The mechanism 
for generating neutrino masses is discussed in Section III for the three different cases and 
compared to the corresponding Oz//3/3 decay amplitude. In Section IV we identify the new 
particle additions generating those operators at tree level after integrating them out. A 
more detailed exposition is presented in the appendix. Whereas explicit (renormalizable) 
models are given for each scenario in Section V. The phenomenological implications for 
LFV processes and LHC searches are reviewed in Section VI. The last section is devoted to 
conclusions. As indicated, technical details are collected in the appendix. 
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II. LOWEST ORDER EFFECTIVE OPERATORS CONTRIBUTING TO DE- 
CAY 



The effects of NP below its characteristic scale A are adequately described by an effective 
theory involving the light fields and preserving the unbroken local symmetries; the smaller 
(compared to A) the characteristic energy scale of the relevant processes, the higher the 
accuracy of such a description. This justifies parameterizing Oi//3/3 decay, whose effective 
scale is some tens of MeV, by a set of effective operators describing NP above the electroweak 
scale. These operators must violate LN, as does Oz//3/3 decay and the observed neutrino 
masses, which we assume to be of Majorana type (we also assume there are no light RH 
neutrinos). 

In this section we will classify the lowest order local operators mediating Ou/SP decay. As 
mentioned in the introduction, this has been addressed previously in the literature, often 
including quark fields [20-23] (see also [24, 26] for recent models realizing operators with 
quark fields). In contrast to these papers, we shall assume that the NP does not couple 
directly to the quark sector, so that all quark interactions are mediated by the electroweak 
gauge bosons. In the appendix we provide the methodology for constructing the effective 
operators to any given order, discussing in detail the dominant ones for each final electron 
chirality. 

In order to write down the basis of effective operators we first have to fix the light field 
content and the symmetries the operators must satisfy. We shall restrict ourselves to the 
SM fields and local symmetries, although allowing for more that one light scalar doublet, 
as this gives a few more alternatives and may simplify explicit realistic realizations, as we 
shall illustrate later on. Operators contributing to Oz//3/3 decay must involve two leptons of 
either chirality, ii or e^, and a number of scalar doublets, (p, to make the product invariant 
under SU{2)l x f/(l)y transformations. Besides, they can have covariant derivative (-D^^) 
insertions, which do not change the field quantum numbers. We will be mainly interested 
in the lowest order operators of each class, and in the heavy particles whose virtual effects 
can generate them (assuming only renormalizable couplings). 

Before proceeding we note, that in some cases the number of different light scalar doublets 
(pi in the theory matters, since for i > 1 qualitatively new operators become possible; this is 
related to the appearance of additional physical scalars. When more than one light doublet 
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is present one must consider all possible scalar flavor assignments to all allowed operator 
structures. For example, in the RR case we have an operator of dimension 7, analogous 
to 0^^\ but with no covariant derivatives: frRe^(0|0j)^. This field product vanishes if 
i = j because the scalar product of two isodoublets is antisymmetric; however, if we 
have two light scalar doublets, the operator does not vanish. Such a theory will contain 
light charged scalars if^; upon spontaneous symmetry breaking, (pjcpj oc H~ + so 
that operator generates the coupling e^e'j^H^ and will contribute to Oz//3/3 decay if the 
physical charged Higgs H~ couples to quarks. Similarly, the number of sets of possible heavy 
excitations that generate the effective operators after being integrated out increases with 
the number of independent light scalar doublets. 

We will proceed as follows. In the classification of independent operators we will for the 
most part assume there is only one SM Higgs doublet. We will do the same for the listing 
of the SM additions (that is, the heavy excitations) generating these operators; in the ap- 
pendix we also discuss the differences when there is more than one light scalar isodoublet. 
In the model generating O^"^"* below (Section V 2) we make use of two light scalar doublets 
for illustration, and simplicity; in that particular case the presence of a discrete symme- 
try guarantees that there is only one lowest-order effective operator involving two leptons, 
one LH and one RH, despite having two doublets (technically this is because a covariant 
derivative and a gamma matrix are both required to match chirality and Lorentz indices). 



1. LNV operators with two LH leptons 

The lowest-order operator of this type is the only invariant dimension 5 (Weinberg) 
operator [32, 33] displayed in Eq. (2). The corresponding term in the Lagrangian is 

where the isodoublet corresponds to the first generation (as will be the case throughout this 
section). In general there are corresponding contributions for all lepton flavors, see Section 
III, which provide Majorana masses to all the SM neutrinos after electroweak symmetry 
breaking. This operator contributes to Oz//3/3 decay by transforming the two neutrinos into 
two LH electrons by the emission of two W gauge bosons (see Fig. 2). 

All LNV operators of dimension 6 involve quarks and violate baryon number by ±1 unit; 
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FIG. 2: Tree-level diagram mediated by O^^-* (light neutrino masses) contributing to 0^(3(3 decay. 

they also conserve B — L and will not contribute to Oz//3/3 decay. There are, however, three 
independent operators of dimension 7 that do contribute to Qvj3(5 decay: 

(see the appendix). One may think that the operators in Eq. (9) generate important 
contributions to LNV processes; this, however, is not the case. The reason is that any 
model generating (9) at tree level necessarily also generates (9^^^ in Eq. (8) at tree level; all 
these operators contribute to all the processes we are interested in, but those from (9) will 
always be suppressed by two additional powers of A and are subdominant. Accordingly, we 
will ignore these operators in the phenomenological analysis. 

The limit on Oz//3/3 decay, for instance, from ^^Ge [40, 41] (see also [11, 12] for recent 
reviews) is usually expressed in terms of the electron-electron element of the neutrino mass 
matrix 

|(m^)ee| < 0.24-0.5 eV. (10) 

This constraint on {my)ee is consistent with the neutrino mass limits from oscillations, 
cosmology and tritium beta decay. When (8) is used to express {mv)ee in terms of Cie"*, the 
following restriction on A is obtained 

^ >10"TeV. (11) 



\cfj\ 

The bound in (10) can be translated into a limit on the amplitude for Oi//3/3 decay at the 
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FIG. 3: Tree-level diagram contributing to decay with one insertion. 



q2 
Peff 



(12) 



parton level which can be estimated as 

1-^0 

where PeS ~ 100 MeV is the neutrino effective momentum obtained from averaging the 
corresponding nuclear matrix element contribution. Thus, from (10) one obtains 

-9 



Q2 l-^Oi'/?/?! — ^ O X iU 



(13) 



•F Pes 

which is the limit on the Oz//3/3 decay amplitude that we will also impose in the other two 



cases. 



2. LNV operators with one LH lepton and one RH charged lepton 

The leading operator of this class, given in Eq. (3), has dimension 7 because it must involve 
two leptons, three scalar doublets to cancel the leptonic hypercharge, and one covariant 
derivative to compensate a chirality flip. Explicitly, 

0£) = {ct^^D^~^)ct^^e^^;U ^ i^v'W;:e^^^v:^ + ■■■ . (14) 

It must be noted, however, that in order to write simple, working models fulfilling all the 
experimental requirements, it may be necessary to impose additional symmetries, and con- 
sequently, some of the external (light) scalar doublet fields may not coincide with the SM 
Higgs doublet, as we mentioned above. This will be the case in the explicit model we will 
work out below. 

The estimate from the Oz//3/3 decay amplitude shown in Fig. 3 is given by 
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FIG. 4: Tree-level diagram contributing to decay with one insertion. 

This translates into a bound on A, which must be > 100 TeV, if we want A^JJj^p 

to satisfy the hmit in (13). These are order of magnitude estimates, but one can also 
use detailed nuclear matrix elements available in the literature [42, 43]. The interaction 
induced by the operator O^'^'^ can be partially expressed as a modification of the standard 
weak interaction, W^e^y^ ((1 — 75) + 77(1 + 75)) u, where u = + is a. Majorana field. 
Then, the strong limit on rj derived using detailed nuclear matrix elements calculations, 
\r]\ < 4.4 X 10"^ (see [42] and [43] where rj was termed Cyl^), reads in our case 

l^l = |^|C^i?l<4.4xlO~^ (16) 
implying a bound which is very close to our estimate 

^ > 106 TeV. (17) 



3. LNV operators with two RH charged leptons 

In this class the leading operator has dimension 9 and it is given in Eq. (4): 

C^ie^ = e^el^{<P^Dy){<P^D,4>) -^-^v'W;;W-^e;^el^ + ■■■. (18) 

As we have argued in the companion paper [35], the 0^(3(3 decay amplitude in Fig. 4 can 
be large, near to its present experimental limit, with 

10 
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LL 


LR 


RR 






10^ 
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TABLE L Natural (with cie^ ~ 1) NP scale limits for the three different lowest order effective 
operators mediating Oi//3/3 decay. 

Requiring that the ^q^'*^^ amphtude satisfies the bound in Eq. (13), one obtains A > 
2|deV^^ TeV, but one can also make use of the detailed nuclear matrix element calcu- 
lations. From O^^^ one obtains the following six-fermion contact interaction inducing Ou(3(3 
decay 

Cou(sp = (uYil - 75)rf) (^7m(1 - Ts)^/) e(l - 75)e'^ , (20) 

where nip denotes the proton mass and in our case 

es • (21j 

This type of interaction has been studied in [44], where the bound |e3| < 1.4 x 10~^ at 
90% C.L. was found ^. This also leads to a limit very close to our estimate 

> 2.7 TeV. (22) 

In Ref. [35] we present a realistic model which accommodates the observed neutrino masses 
and a large Oz//3/3 decay observable in the next round of experiments, proving explicitly 
the consistency of the previous estimates. In particular, neutrino masses are naturally 
suppressed even with NP at the TeV scale. 

These scale estimates are summarized in the Table I. Note, however, that in actual 
models the operator coefficients d? are not in general ~ 1 and the A estimates may vary. 
For the LL case there are models with NP at A ~ 1 TeV for sufficiently suppressed couplings 
C*-^-* ~ 10^^^ (see, for instance, Refs. [45-49]). We can also have different scales within the 
model, as in the LR model below, where the new leptons can have masses below a TeV and 
be observable at the LHC Similarly, in the RR case in Ref. [35] the new scalar masses can 
range from few hundreds of GeV to tens of TeV. A more detailed discussion can be found 
in the companion paper. 

^ There is a misprint in Ref. [44]. We thank the authors of this reference for providing us with the correct 
limit on 63. 
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III. MAJORANA NEUTRINO MASSES GENERATED BY THE LNV OPERA- 
TORS INDUCING 01//3/3 DECAY 



Once the effective theory generates one of the LNV operators that produce Oz//3/3 decay, 
neutrinos will get a mass at some (loop) order, even if there is no other independent source 
of neutrino masses. The three operators not only stand for NP at quite different mass scales 
but result in different neutrino mass structures. Specific models may of course include a 
combination of these effects, so that the final structure may be quite involved. 



1. LL operator 

After electroweak symmetry breaking O^^^ generates a well-known and much studied 
contribution to the neutrino masses (summation on repeated indices must be understood 
when applicable through the manuscript): 



^(5) ^(5) „2r'(^) 

jy ab jy 



aL<p(P'hL 



A 



(23) 



where a,b = e, fj,,T are flavor indices. Hence, 



[m^ ab 



A 



(24) 



In this case Oz//5/3 decay proceeds through the diagram in Fig. 2, and it is proportional to 



(25) 



where U is the mixing matrix diagonalizing the neutrino mass matrix written in the current 
eigenstate basis with well-defined charged lepton flavor [4, 5]: 



2t;2C(5) 



t 



U 



nil 



\ 



m2 



and 



U 



C13C12 



Cl3'5l2 



si^e 



-C23-512 - S23Si3Ci2e*'^ C23C12 - S23-5l3Sl2e*^ S23C13 



S23S12 - C23Si3Ci2e' 



-S23C12 - C23Sl3Si2e*'' C23C13 

12 



/2 



(26) 



3102/2 



(27) 



where Sij = sin 6ij and Cij = cos%. Thus A is required to be ~ 6 x 10^^ TeV to reproduce 
the observed neutrino masses {m^)rT ~ 0.1 eV for Cj^f^ ~ 1, which is of the same order as 
the hmit derived from Ou/3f3 decay (see Table I). 



2. LR operator 

Restoring flavor indices the LR contribution to the effective Lagrangian becomes 

-^{(P^DyWe^lJbL) . (28) 

In this case the neutrino masses are generated by radiative corrections (left graph in Fig. 1). 
Here it is important to differentiate between the calculable (logarithmic) contributions to 
the masses derived from O^'^^ and the estimates obtained by matching. The first, obtained in 
standard effective field-theory fashion using dimensional regularization and a renormalizable 
gauge, scale like 



{6m 



y3 



u lab ^ 



(m^dJJ + m.cil^) log . (29) 



167r2A3 

Whereas the estimates from matching, which give the dominant contribution to neutrino 
masses, are obtained using dimensional analysis and are of the form 

V 



[m^ ab 



(m„Cl? + m.cl'j) , (30) 



levr^A 

as we will derive for the specific models discussed below. Note, however, that even for this 
last estimate it is important to use a renormalizable gauge as in the unitary gauge spurious 
positive powers of A may appear. That the final result must be proportional to 1/A follows 
from (i) the fact that the dimension 5 Weinberg operator is the only one describing neutrino 
Majorana masses (in the absence of uji); and (ii) that we assume that the NP is decoupling. 
It is also important to note that, as mentioned in the introduction, specific models may have 
more than one NP scale so that A represents an effective scale and may not correspond to 
the mass of any specific particle. For example, if the theory has two scales M > M' one 
may have A = M'^/M'; we will present an example of this situation in Section V 2. 

The expression (30) coincides with the results obtained in specific models such as the 
one worked-out in Section V 2 (Eq. (41)) with the appropriate identification of A. In this 
case {m^)rT ~ 0.1 eV with ciV ~ 1 implies (see Eq. (6)) A ~ 4 x 10^ TeV, that should 
be compared to the limit obtained from Ou/^P decay in Table I, which is several orders of 
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magnitude smaller. This implies that in realistic models the coefficients C*^^^ must be much 
less than 1 to allow for a A of the order of 100 TeV, in agreement with the Oz//3/3 decay 
estimate, if this is to be observed in the next generation of experiments. This is what 
happens in the explicit model we will work out below. 

3. RR operator 

Finally, the LNV operator (9^^^ generates Majorana masses to neutrinos at two loops, in 
this case also suppressed but by two loop factors and two charged lepton mass insertions 
(see the right panel of Fig. 1). As for the LR case we distinguish between the calculable 
(logarithmic) contributions, 

^ (i6j2)2A5 ^°^a?^fc log , (31) 
and the estimates obtained from matching 

Explicit calculations in specific models (using renormalizable gauges) [35] reproduce (32) up 
to a proportionality factor of order 1 that depends on the various masses in the loop, and 
when A is identified with an appropriate combination of heavy scales. For [m^)rT ~ 0.1 eV 
and Ci^T ~ 1, (32) implies A ~ 1.3 x 10^ TeV, several orders of magnitude larger than the 
0iy(3(3 decay estimate in Table I. Although in realistic models the coefficients C-^-* and A are 
typically smaller [35]. 

IV. HEAVY PARTICLE ADDITIONS GENERATING THE LOWEST ORDER 
LNV OPERATORS AT TREE LEVEL 

In this section we will work out the combinations of heavy particles (scalars, fermions and 
gauge bosons) that must be present in any extension of the SM if it is to generate one of the 
operators O^^'"^'^^ at tree level; thus producing the largest possible rates for 0iy(3(3 decay. We 
will assume that the underlying theory is weakly coupled and contains only renormalizable 
vertices ^; we also assume the NP respects all the gauge symmetries of the SM. In listing the 

^ Non-renormalizable vertices are presuraably suppressed by inverse powers of a yet higher scale Ahigh that 
we assume much larger than A. 
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FIG. 5: Topologies generating 0^^\ These graphs together with those with one or two W boson 
attachments generate 0^'^'^''^'^^''^'^^^K 

(Y) (Y) (Y) 

heavy particles we denote by X) , and $} a heavy vector, fermion or scalar with 

isospin I and hypercharge Y, respectively. When the heavy particles can be either a heavy 

(Y) 

vector or heavy scalar with the same isospin and hypercharge, we use 5) to denote both 
possibilities. 



1. LL additions 



In the case of LNV operators with two LH leptons the heavy excitations that can generate 
the operators of dimension 5 and 7 in Eqs. (8) and (9), respectively, are 

v,(0). (33) 

(If the model has two or more light scalar isodoublets, a heavy scalar singlet $q^^ can also 
generate these operators [50].) This means that, as mentioned previously, if the underlying 
theory generates (9), it will also generate (8) with the same heavy scale. In Fig. 5 we depict 
the diagram topologies resulting in those operators upon integration of the heavy particles 
flowing through the internal lines. 



2. LR additions 



The sets of heavy excitations that can generate the dimension 7 LNV operator with one 
LH lepton and one RH electron in Eq. (14) are 

ry{3/2) . r,T/(V2) . JxT/(l/2) y (0). . r -. (0) ^(l)-,. r,T,(0) fnA\ 
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In this case one must always exchange two heavy particles. It is worth noting that in the last 
two possibilities the heavy fermions \l/g''J necessarily have couplings that would also generate 
the operators in Eqs. (8) and (9) at tree level. In contrast, models containing $g^] may or 
may not generate them, depending on whether the heavy scalars couple to £l x (which 
may be forbidden by the symmetries of the underlying theory). In the appendix we provide 
the diagram topologies relevant to each case. 

One can construct many models choosing from the above NP matter contents. However, 
in order to make the model phenomenologically viable we must often enlarge these minimal 
sets. This is because (34) are fixed only by the requirement that they generate O^'^'^ at 
tree level, which does not insure the preservation of extra symmetries that are sometimes 
necessary, for instance, to avoid too large LFV rates, or to forbid tree-level neutrino masses, 
or to reproduce the observed lepton spectrum. In Section V we discuss a realistic, simple 
example for the case But, as we will argue, for the model to be reahstic, it 

must include at least two fermion doublets ^[/2^^ besides one scalar triplet \ and, in 
addition, a second light scalar doublet 0'. 

3. RR additions 

Finally, the sets of heavy excitations that can generate at tree level the dimension 9 LNV 
operator with two RH electrons in Eq. (18) are 

{*;'//',<'}; {*S°',b|"}; {*r,B»'}, (35) 

Despite the presence of ^l^l and in some of these options, these heavy particles need 
not have the same vertices as the ones leading to (8). If they do, O^^^ would have only 
subdominant effects; but this is in general not the case. In the appendix we also provide the 
diagram topologies relevant to each of these cases. 

There are many models that can be constructed containing the above particle content. 
For example, the case where the scalar sector of the SM is extended by adding a doubly- 
charged isosinglet $0^^ and an isotriplet of unit hypercharge \ was considered in detail 
in the companion paper [35]. (See also [51, 52].) As in other cases, additional heavy fields 
may be required in order to make the model realistic. 
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V. SIMPLE EXAMPLES OF FUNDAMENTAL THEORIES WITH A LARGE Oi//3/3 
DECAY RATE AND REALISTIC NEUTRINO MASSES 



There are many specific models fulfilling the three generic scenarios discussed above 
(see Eqs. (2), (3) and (4)). Here we list a few for illustration purposes; we do not aim at 
reviewing all models that have been considered in the literature; concentrating instead on 
specific examples that exhibit the salient features discussed previously. In particular neutrino 
masses take the form in Eqs. (5), (6) and (7), and in general the coefficients C^^^ <^ 1 for 
consistency in realistic models. 

1. LL models 

SM extensions that at tree level generate the Weinberg operator O^^^ at low energies, 
have been extensively studied in the literature (for type I see-saw see [27-30], for see-saw 
type II see [53-58], and for see-saw type III see [59, 60]). Specific models for any of the 
three possible see-saw scenarios ^ (see Eq. (33)) require very heavy mediators, with masses 
~ 10^^ GeV or very small couplings, and were discussed some time ago (for recent reviews 
see [2, 3, 61, 62] ). The most popular and simplest case, and also the pioneering one, results 
from the addition of heavy RH neutrinos '^'^^ = ur with the renormalizable Lagrangian 

-C""^ = iT^$l'aR - {^MaU^R VaR - yabT^4>^ ^bL + h.C.} , (36) 

including the kinetic terms and the Yukawa couplings Uab- There must be at least two 
heavy neutrinos VaR (with masses Ma) to guarantee that the light neutrino mass matrix is 
at least of rank 2, required in order to account for the three non-degenerate light neutrinos. 
Explicitly, 

< 10"" TeV"^ , (37) 

in agreement with Eq. (11) and Table I. The alternative case with the see-saw messengers 
near the TeV scale (and \y\ < 10~^) has become more popular with the launch of the LHC, 
see [48, 63], and references therein. 

In this scenario Oz//3/3 decay follows from the exchange of the light Majorana neutrinos; the 
resulting amplitude is proportional to the effective electron-neutrino mass (mi,)ee in Eq. (25), 

^ Assuming only one light SM Higgs doublet. 
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see Fig. 2; hence, any Oz//3/3 decay rate within the present experimental precision can be 
accommodated. Indeed, a global fit to neutrino oscillation data gives (see, for instance, [6]) 
Amii =ml-ml = {7 x IQ-^ eV^, Am^^ = ml - ml = (2.50;°;?^) x IQ-^ eV^, 
s?2 = 0.312j:°:°}^, = 0.52+°:°^, = 0.013+|];[{|J^. Replacing these values into Eq. (25) one 
can obtain at Icr any Oz//3/3 decay rate compatible with present experimental limits for the 
normal hierarchy, although it is bounded from below for the inverse one [1, 64, 65]. 

2. LR models 

To our knowledge, no realization of this second scenario has been spelled out in the 
literature. One can construct many models with Oz//3/3 decay into two electrons of opposite 
chirality mediated by O^''^ by choosing among the matter contents in Eq. (34); however, 
when constructing a realistic model we must in general enlarge these minimal sets. 

We start from the set {^'^y^^^^^i^} contained in Eq. (34); to simplify the notation we 
define = a scalar isotriplet of hypercharge 1, and ^I^^/2^'' = = + L^, a lepton 
isodoublet of hypercharge 1/2 (in terms of its LH and RH components); a simple way to 
insure the decoupling of the heavy physics is to assume, as we do, that the heavy fermions 
are vector like. This particle content is sufficient to generate O^"^^ at tree level, and it is 
not hard to convince oneself (see the appendix) that the relevant graphs must involve the 
couplings eR(f)L, iiLx and However, such a model also allows the coupling Il^lX and 

will then generate O^^^ at tree level through the standard type-II see-saw diagram (on the 
right of Fig. 5), a possibility we wish to disallow. In order to do this we impose a discrete 
Z2 symmetry under which x and L are odd and £l is even, so O^^^ does not appear at tree 
level; unfortunately, this symmetry also forbids the vertex . In order to overcome this 

difficulty we assume the presence of two light scalar doublets 0, 0' which are, respectively, 
even and odd under Z2. The allowed vertices are then 

Lxh, L^'cR, 0^x0', (38) 

from which O^"^^ is generated through diagrams such as the one in Fig. 6, while O^^^ appears 
only at one loop. 

As we will explain below, in order to accommodate a generic neutrino mass matrix and 
also allow for flavor symmetries treating the three families on the same footing, we will 
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FIG. 6: Tree-level diagram contributing to O*-^) in the model proposed. 
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TABLE II: Quantum number assignment for the extra fields. 

assume the presence of 3 heavy vector-hke fermion doublets L^, a = 1,2,3. The complete 
list of new fields is given in Table II; the Lagrangian will include renormalizable couplings 
preserving these symmetries, noting that the SM fields transform trivially under Z2. The 
second scalar doublet 0' could be identified with the isodoublet giving mass to the up quark 
sector in two-doublet models [66-71] (for reviews see for instance [72, 73]) if we require the 
RH up-quark singlets to be odd under Z2 (the quark isodoublets are even). We will prefer 
not to do so, because we will find it convenient to be able to assign a small VEV (0'). We 
also note that in this model LN is explicitly broken by (renormalizable) terms in the scalar 
potential, in particular by (0"'"0')^. 

The scalar potential can be easily arranged to insure a minimum where (0) ^ {</)') , (x) 7^ 
0, with (x) — — /U* (0') (0) /m'^, /i the trilinear 0^X0' coupling and the isotriplet mass 
(in order to satisfy the limit from electroweak precision data [1, 74] we require (x) ^ 2 
GeV; see [35] for a similar analysis and [75] for a recent study, in the framework of the 
type II see-saw, including one- loop radiative corrections from the scalars of the model). We 
assume negative mass terms for and 0' to trigger the corresponding VEVs, whereas x g^ts 
a VEV through its mixing with the scalar isodoublets. Otherwise, dimensional couplings 
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in the potential are typically of electroweak order, except for new scalar masses that may 
be larger. Dimensionless ones stay perturbative, in general ranging within an uem ~ 10~^ 
factor. 

We now discuss briefly the viability of the model, concentrating on the effects of the 
masses and mixings of the new fermions, the induced LFV effects, and the implications for 
both the LHC, and the light neutrino masses. The heavy lepton Lagrangian reads 



£^ = La{i0- Ma)La + {y^feL^i^'cfej? + y'^.L^LXhi + h.c.} , (39) 

where we assumed that the heavy mass matrix is diagonal without loss of generality. Once 
(p' and X acquire VEVs the light e^, ii leptons mix with the L^, such mixings and the 
corresponding phenomenology of heavy vector-like lepton doublets were analysed long ago 
[76-78], and more recently within the context of Little Higgs models [79-82], and Extra 
Dimensional theories [83-85] (for a review see [86, 87]; for updated limits see [88]). 

The low-energy effects of such mixings will be proportional to (0') /Ma or y^^, (x) /Ma 
and can be made as small as experimentally required by increasing the heavy masses Ma, 
reducing the couplings y'^^, or the VEVs (0') , (x)- LEV effects can be further suppressed 
by assuming that the light charged leptons, which get their masses through the SM Higgs 
mechanism, are aligned along the heavy flavors. This corresponds to taking diagonal, 
which may be natural in a larger model. The LHC reach for the scalar triplet was reviewed 
in the companion paper and is updated in next section. The production of the La at the 
LHC has been also studied previously [89] and the general conclusion is that they will be 
detected provided their masses are below 850 GeV for a center of mass (CM) energy of 
14 TeV and an integrated luminosity of 100 fb~^ [90] (the LHC reach reduces to 350 GeV 
for heavy leptons mainly decaying into taus [91]). 

(7) 

Given the couplings of the model one can evaluate C^^ by using the diagram in Eig. 6, 

A3 ~ m2 M2 ' ^^^^ 

where all masses in the L^. and x multiplets are taken equal. 

Light neutrinos are massless at tree level, but they get a mass at one loop through 
diagrams like that in Eig. 7. The full calculation gives 

^'"'^^'^^ ^ 32^^ (^acc + ^bv^y::) log ^ , (41) 

20 





\ 

\ 


t 

t 


1 
1 


\ 
1 









FIG. 7: Leading one-loop contribution in the Feynman gauge to neutrino masses in the fundamental 
theory. 

where v' = {</)') and we have assumed that all other masses are much smaller than Mc and 
m^. Thus, with only one heavy lepton doublet the neutrino mass matrix has at most rank 2. 
With two heavy lepton doublets all three light neutrinos can be massive, but with three it 
is also straightforward to impose a flavor symmetry that forbids any potentially large LEV. 

The model may face a domain wall problem if the spontaneously broken Z2 symmetry is 
exact. This can be obviated adding a softly breaking term to the scalar potential. In 
a future publication we will provide a detailed analysis of this model including the neutrino 
mass calculation, as well as the quantitative discussion of the parameter space allowed by 
present experimental bounds, and the predictions for the different observables. 

3. RR models 

The final case, where Oi//3/3 decay involves two RH electrons through the operator 0^^\ 
allows for many tree- level realizations, typically with a heavy sector near the electroweak 
scale as argued in Section IV and in the appendix. A realistic, simple model of this scenario 
with the neutrino masses generated at two loops, and predicting a non-zero third mixing 
angle sin^ ^13 > 0.008 for a large Oz//3/3 decay rate, was studied thoroughly in the companion 
paper [35] and corresponds to the particular SM addition {^'^\ ^^i^} in Eq- (35). The model 
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contains in addition a real scalar singlet a whose presence allows a Z2 discrete symmetry that 
protects the neutrinos from acquiring a tree-level mass. This, in turn, makes the neutrino 
masses calculable. 

VI. OTHER PHENOMENOLOGICAL IMPLICATIONS 

Although it is beyond the scope of this paper to make a detailed discussion of the exper- 
imental constraints on the different scenarios and models, we will briefly comment on the 
main phenomenological implications for ongoing and forthcoming experiments. The present 
study answers the general question of how large can the Oi//3/3 amplitude be, requiring that 
the only light fields are those of the SM, and assuming that the NP responsible for the 
effect is weakly coupled and has a characteristic scale above the electroweak scale. In par- 
ticular this means then that neutrino masses are Majorana; we also assume, as it is widely 
believed, that the observed pattern of neutrino oscillations indicates the presence of three 
massive neutrinos [1]. This strongly influences our analysis because we must then explain 
the neutrino spectrum, assuming that there is no other larger source of LNV beyond the 
one mediating Oz//3/3 decay. 

Thus, within this framework, we have to verify on a case by case basis whether: 

• The ^1/(3(5 decay rate can be fast enough to be observable at the next round of exper- 
iments. 

• The neutrino masses are correctly predicted. 

• LFV and universality limits are within experimental bounds. 

• The new heavy particles satisfy the collider exclusion limits. 

We have already addressed the first two points in general and in some detail; we consider 
the remaining two in the following subsections. 

A. LFV processes in SM extensions with a sizeable Oi^/3/3 decay rate 

A relatively large Oz//3/3 decay rate requires a not too heavy NP (except for the case 
where O^^^ is generated at tree level). This, in turn, implies that the effects of the new 



22 



particles may be detected in highly suppressed processes, like in flavor-changing leptonic 
or Z decays [1]; the corresponding constraints are model dependent, however. We can 
compare, for illustration, the LR and RR models in the previous section and in Ref. [35], 
respectively. In the latter the heavy sector only involves new scalars, with little effect on 
lepton universality precision tests, for instance. The main universal constraint derived from 
them being the upper bound on the isotriplet VEV, (x) ^ 2 GeV [1, 74]. This limit is easy 
to satisfy although it is rather restrictive in this speciflc model due to the small number of 
new free parameters available, and the desirability of accommodating a large Oi//3/3 decay 
rate. For the same reason LFV constraints are quite demanding because the corresponding 
leptonic decays are proportional to neutrino masses, which have a large misalignment from 
the charged lepton current eigenstates, described by the Pontecorvo-Maki-Nakag aw a- S akat a 
mixing matrix [4, 5]. This translates into lower bounds on the heavy scalar masses with a 
preference for values above a TeV [35]. 

The LR model presented in this manuscript is somewhat different, for the new scalars 
have no direct coupling to SM lepton pairs, although they are induced by fermion mixing. 
In this case the most stringent restrictions concern the new heavy leptons, with masses M^, 
especially on their mixing with the light ones; such indirect constraints on extra vector-like 
leptons have been thoroughly studied in the literature [76-88]. The vector-like character 
that allows their decoupling without breaking the SM, also flxes their mixing behavior 
and their low energy phenomenology. The corrections to SM vertices are suppressed by 
at least 2 powers of a small ratio 0{m^/Ma)i where is typically a light lepton mass. 
Hence, the decay rates, which are vanishingly small within the SM, are suppressed by at 
least 4 powers of these ratios; moreover, they also vanish when (0') — )■ 0. Even more, 
LFV processes can be also canceled by aligning the heavy and light lepton flavors. In 
summary, many different small factors can conspire to make negligible the effects of the heavy 
fermion, in particular, their tree-level effects are naturally small enough to accommodate the 
experimental constraints. The main restrictions on this type of models result from one-loop 
contributions exchanging heavy leptons and bosons; the most restrictive processes being 
those including the muon to electron transition. In fact, one can, to a large extent, apply 
the conclusions from related analyses for the Littlest Higgs model with T-parity [79-82]; the 
general conclusion is that the heavy flavors must be aligned with the light charged leptons 
with a precision better than 1 — 10 % for heavy masses of O(TeV). 
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In summary, LFV provide stringent restrictions on these models but they can be satisfied 
within relatively large regions of parameter space. What would be more interesting, these 
models could also explain a departure from the SM predictions if found in the ongoing and 
forthcoming searches for LFV [1, 92] (see [93] for a review). 

B. Collider searches for particles with LNV interactions 

The specific models above also have different collider signatures: the LR model contains 
heavy fermions and scalars, while the RR model contains only heavy scalars; all of which 
can be searched for at the LHC. Verification of either model would involve not only the 
discovery of the corresponding heavy particles, but also a demonstration of the presence of 
LNV interactions, (the possibility of observing LNV events at LHC was emphasized quite 
some time ago [94]). Although not all decays of the new particles produce LNV signals, 
in some cases they may be dominant; however, in general one has to search for the new 
resonances in the most sensitive channels, and only afterwards address the possibility of 
observing LNV events. Typically, these will be difficult to observe because LNV is usually a 
small effect since the corresponding amplitudes involve several (small) couplings all of which 
must be present in order for LN to be broken. Hence, in general the dominant production 
mechanisms are standard and LN conserving, otherwise they are small (though LNV decays 
can be slow, they typically still can occur within the detector [95]). 

Following the discussion in Ref. [35] we will first comment on the detection of the new 
scalars in the RR model. Doubly-charged scalars have fixed couplings to photons, so that 
their production cross section is known; their decay into leptons (if allowed) gives a very clean 
signal, which is particularly important at hadronic machines; although it is not LNV by itself. 
Therefore, if doubly-charged scalars are light enough, they are quite suited for detection at 
colliders. Generally, this type of scalars forms part of a weak triplet, and usually also acts as 
see-saw messenger of type II [53-58, 96] (see also [97]). These triplets are then theoretically 
well-motivated, especially when considering L-R symmetric models, and simulations of their 
production at future colliders can be found in the literature [98-102] (see also [103, 104] for 
recent studies; and for model-independent ones [105, 106]). The general conclusion is that 
the LHC discovery limit can reach masses over 600 GeV (for a CM energy of 14 TeV and 
an integrated luminosity of 30 fb~^) [103, 107]; although the actual limits may be larger 
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given the outstanding LHC performance for a CM energy of 7 TeV [108]. (See for a recent 
review [109].) First results from CMS [110] and ATLAS [111] have been presented at this 
last CM energy with an integrated luminosity of 0.89 fb~^ for any di-lepton scalar final state 
and of 1.6 fb^^ for di-muon final states, respectively. As nothing is seen, a lower bound on 
the doubly-charged scalar mass is obtained: of about 250 GeV if the main decay channels 
contain r leptons and to about 300 GeV if they contain only electrons or muons [110], 
reaching 375 GeV if they only couple to muons [111]. (Present Fermilab Tevatron Collider 
limits are less stringent [112, 113].) In our case, however, the doubly-charged triplet, 
does not directly couple to fermions; while the other doubly-charged scalar in the model, 
a singlet k^^, does not couple to W pairs. However, they mix; and both of them can be 
produced at LHC via the Drell-Yan mechanism {qq — )■ 7*,Z* — )■ x~^^X ). Since 

this is the main production process assumed by both LHC Collaborations, the former limits 
apply directly to the singlet decaying dominantly to lepton pairs for a small mixing: > 
300 GeV. Limits on the triplet mass will be more difficult to derive, for the dominant process 
qq — )■ 7*,Z* — )■ x~^^X ~^ W^W^W~W~ is more complicated to study, due to its large 
backgrounds and the inherent difficulty of reconstructing several leptonic W decays [114]. 
The only viable LNV decay channels qq — ?■ 7*, Z* — )■ x^^X 1 i^^^f^ e^e^W^W^ are 
suppressed by the small mixing between and k^^. 

The LR model discussed above contains, besides the scalar isotriplet of unit hypercharge 
X, a second scalar isodoublet and several vector-like lepton doublets. This scalar triplet 
couples to a light and a heavy lepton, and the latter decays into a light lepton and a W, or a 
Z, or a Higgs boson; this results in a four-fermion decay. Similarly, the extra scalar isodoublet 
decays into four fermions. Then, if also pair produced via the Drell-Yan mechanism, final 
states will have at least eight fermions. The signal may be striking due to the large number 
of charged leptons, but there are many open channels and may be not easy to resolve the 
different samples. These final states are different from those of the RR model, which will 
eventually allow to discriminate between both theories. 

The heavy vector-like lepton doublets are also mainly produced in pairs. They violate 
the Glashow-Iliopoulos-Maiani mechanism [115], and can decay through a flavor changing 
neutral current into a light lepton and a Z or Higgs boson [89]; as for sequential fermions they 
can also decay into a lepton and a W boson through the usual charged current interaction. 
In these cases the flnal states have at least six fermions, so that the heavy leptons will 
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be relatively easy to find if light enough: Ma ^ TeV for a CM energy of 14 TeV and 
an integrated luminosity of 100 fb~^ [89-91]. The dominant decays of these (quasi) Dirac 
fermions are LN conserving to a large extent. 

VII. CONCLUSIONS 

To date there is no direct evidence that LN is not a symmetry of nature [1], at least at 
the energies and experimental sensitivities currently available. There are, however, scenarios 
that suggest LN may be violated at higher scale. For example, the baryon asymmetry of 
the universe (see for example [116]) can be explained through leptogenesis [117], yet no 
experimental indication of such a mechanism has been observed. Light neutrino masses 
and mixing angles provide the most favored explanation for neutrino oscillations, but the 
experiments do not distinguish between Dirac and Majorana masses, and thus, they do 
not require LNV; although, in many models light neutrinos are assumed to be Majorana 
fermions, implying that LNV is also assumed. 

In contrast Oz//3/3 decay and appropriate signals at LHC will be sensitive to LNV effects; 
and in fact, these are the only way known to experimentally establish the presence of LNV. 
Hence the relevance of these two types of experiments. A new generation of Oz//3/3 decay 
experiments are underway (see [11, 12] for recent reviews), where a positive signal would 
provide conclusive evidence of LNV, and would open a new era of experimental searches and 
theoretical studies aiming at isolating the type of NP that mediates such a process. 

The effective Lagrangian approach allows to address this theoretical question with gener- 
ality and has been the subject of the above discussion. As described in the text, and detailed 
in the appendix, we have constructed all gauge-invariant effective operators of dimension 
< 9, violating LN by 2 units and involving two SM leptons (but no quarks) and any number 
of Higgs isodoublets and covariant derivatives (the appendix also lists LNV operators with 
more than 2 leptons). For each of the three possible scenarios; (i) with two LH leptons (LL), 
(a) one LH lepton and one RH electron (LR), and (in) two RH electrons (RR), there is only 
one lowest order effective operator, 0^^\0^'^'> and 0^'^\ respectively (see Eqs. (2), (3) and 
(4)). They describe the largest possible contribution to Oz//3/3 decay for each final electron 
chirality assignment. We have also identified the possible new particle additions that can 
generate these operators at tree level (assuming that the full theory is renormalizable). 
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In general, models of neutrino masses can have an origin different from that of OuP/S decay; 
irrespective of that, once the NP generates any of the LNV operators, the theory will generate 
neutrino masses at some loop order: at tree level, at one loop or at two loops, depending on 
whether the operator is 0^^\ O^"^^ or 0^^\ respectively. For O^"^'^^ the corresponding masses 
are relatively suppressed by loop and light mass (generated by chirality flips) factors. In 
the case of C(^) there are no such suppressions and the same parameter giving the effective 
neutrino electron mass {m^)ee enters in the 0^(3(3 decay amplitude (see Eq. (12)). Thus, 
both pieces of data only constrain the \C^^^\/A ratio. In contrast in the other two cases, 
and since the neutrino mass scale is fixed to be around 0.1 eV, the observation of Oi^/3/3 
decay would allow to estimate the scale A of NP and the corresponding effective operator 
coefficient C for each scenario for natural theories (and perturbative couplings). Thus, 
leaving to experiments searching for LFV and for collider signatures of LNV mediators. 
To provide an existence proof we have also constructed a simple, realistic model for each 
scenario, explicitly calculating the neutrino masses; as expected from dimensional arguments 
these masses are always proportional to 1/A (see Eqs. (24), (30), (32), (41) and Ref. [35]), 
as argued in the introduction: Eqs. (5), (6) and (7). Then, whereas C^^^ and its associated 
A can vary in between eleven orders of magnitude, at least CtV (CtV) must be at the 
per million (mille) level and the corresponding A of the 10 (few) TeV order. Besides, any 
spectrum of neutrino masses within present experimental limits can be accommodated in 
the first two cases, but in the RR model the neutrino hierarchy must be normal and the 
third mixing angle sin^ ^13 > 0.008, in agreement with recent observations [118-122]. On 
the other hand, once an explicit model is at hand, we can check if it does satisfy the present 
LFV constraints and bounds from large colliders; this is the case for the models discussed. 
Although the Higgs searches now underway at LHC will stringently restrict these models 
with an extended scalar sector near the electroweak scale. 
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Appendix: Effective Lagrangian description of NP 

In this appendix we provide for completeness a brief summary of the effective Lagrangian 
approach and describe the procedures we followed in constructing the operators discussed 
in the main text. 

Effective theories are useful for situations where there is a scale gap: some type of heavy 
physics effects contribute only virtually since the available energies are well below the scale 
of these interactions. Technically one differentiates between the case where the underlying 
physics decouples [34] and when it does not. For the first case the low-energy effective 
theory is obtained by a formal expansion in inverse powers of the heavy scale, the existence 
of which is guaranteed by the decoupling theorem [34] (see also [123]). For non-decoupling 
heavy physics the effective theory is obtained as a derivative expansion [15]. 

Here we shall assume that the heavy physics is both decoupling and weakly coupled, so 
that a perturbative expansion is appropriate; the characteristic scale of these new inter- 
actions will be denoted by A. The general parameterization of NP effects using effective 
interactions is valid at energies below A. The procedure is straightforward: one constructs 
all Lorentz-invariant operators involving the light fields and their derivatives and respecting 
the low-energy local symmetries (here, the SM gauge symmetries); the effective Lagrangian 
is then the linear combination of all such operators, where the A-dependent coefficients 
parametrize all possible (weakly-coupled and decoupling) types of heavy physics. In our 
case the NP will violate LN but will not couple to quarks. If the theory underlying the SM 
were known then one could derive the low-energy effective theory coefficients in terms of 
the parameters of the model. It may then happen that some operators will be absent or 
may appear with suppressed coefficients due to some as yet unknown symmetry. Not know- 
ing the correct SM extension, the effective Lagrangian coefficients are treated as unknowns 
susceptible to experimental determination. 

The coefficient of an operator of (canonical) dimension n is proportional to A^~" so that 
the larger the dimension of the operator the smaller its effect; given a finite experimental 
precision this implies that operators with n sufficiently large can be ignored. In addition, 
operators that are generated by heavy particle loops have coefficients that receive a typical 
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loop suppression factor ~ l/(47r)^. It is important to note that any operator that respects 
the local symmetries of the SM will be generated by the NP at some loop level. Whether 
this happens at tree level depends on the operator and the details of the theory underlying 
the SM; it is a simple exercise to determine the operators and types of NP that have this 
property. It is also easy to device types of NP for which all tree-level-generated operators 
(TLGOs) are absent ''. Dominating effects are then associated with the lowest-dimension 
TLGOs contributing to the process at hand. In this it is important to note that whether 
an operator is generated at tree level or not depends on the details of the heavy physics; 
there are, however, operators that are necessarily loop generated by all modalities of heavy 
physics [125]; we call these loop-generated operators (LGOs). It is also worth keeping in 
mind that in most cases the effects from LGOs will compete with those generated radiatively 
by the SM and are often subdominant; exceptions occur when the SM effects are absent due 
to some accidental SM symmetry, such as custodial symmetry or LN. Except for these 
cases LGO effects lie beyond the experimental sensitivity of current experiments ^. There 
is one additional observation that can be used to simplify the effective Lagrangian: if two 
operators O and O' are such that the combination O — O' vanishes on-shell (that is, when 
the classical equations of motion are imposed), then the S-matrix depends only on the sum 
of the corresponding operator coefficients [126, 127], so that one of the operators can be 
omitted from the effective Lagrangian parameterization ^. 

As emphasized previously, the main feature of the processes we will be interested in is that 
they exhibit LNV. The operators of interest have dimensions ^ 5; many of them have been 
enumerated in earlier publications, including operators involving quarks [20-23], though 
such catalogs are not exhaustive. We do not pretend to provide a complete list of operators 
but concentrate instead on those low-dimensional TLGOs with two leptons and any number 
of bosons and derivatives that can be probed experimentally and provide leading effects for 
wide classes of interesting heavy physics. In the following we will assume that the low-energy 

^ For example, if there is a discrete symmetry under which all the SM particles are singlets but none of the 
new particles are [124]. 

^ This means that radiative effects generated by the heavy excitations are too small to be observed, it does 
not preclude direct observation of new particles (provided the energy available is high enough). In these 
cases the effective Lagrangian approach is, however, inapplicable. 

^ Note in particular that this result implies that the experimental sensitivity to O and O' is the same: one 
cannot replace an operator C by an equivalent one O and find weaker limits on O than on C. 
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excitations are those of the SM, with a scalar sector containing rid > 1 doublets (pi, i = 
1, . . . ,nrf of hypercharge 1/2. Although generally we assumed only one (Higgs) isodoublet 
when classifying the possible higher-order operators, it is worth considering the case > 1 
because this allows the presence of additional operators that are absent in the one-doublet 
case for symmetry reasons; we provide examples of this in the tables below. Moreover, as we 
illustrated in Section V 2, the introduction of extra scalar doublets may allow for simpler, 
phenomenologically viable, fundamental theories. In contrast, no additional light fermions 
are assumed to exist, in particular, any RH gauge-singlet fermions (such as RH neutrinos) 
are assumed to be heavy In presenting our expressions we will not display family indices, 
though in general these are present. 

In order to simplify the notation through this appendix the LH lepton isodoublets are 
denoted by i and RH lepton isosinglets by e; here we will consider operators involving only 
leptons, those involving quarks will be discussed in a future publication. We now provide 
the list of LNV operators of dimension < 9. In doing so we merely provide field content, 
with the understanding that all possible gauge and Lorentz contractions are to be counted 
(so that each of the entries represents, in general, more than one operator). Also to make 
the notation simpler and clarify the physical effects of the operators we find it convenient 
to introduce the following composite operators 

= (l>rL , = ^iD^i^ , J^, = LDn, , = (1>\D^~^, , (42) 

where, as before, 

'''^1, i=eCf, 0. = e0* (43) 

with e = i(T2- Lower-case indices a, 6, etc. are family indices. For the case of a single leptonic 
family and one scalar isodoublet with (0) = (0,f)^ the above operators become 



N = -vul + ■ ■ ■ , 

^9 



+ 



10 



Though it will not be considered in this appendix, it is straightforward to extend the light scalar sector 
by adding a number of light scalar gauge- singlets. 
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= -i^v^W-" + ■■■ , (44) 
v2 

where the elhpsis denote terms involving the physical scalars and we defined ul = CVl^] 
numerically ~ 174 GeV. 

Given an operator O it is straightforward to determine whether it is loop generated or 
whether there are models where it appears at tree level [125]. For the second case one can 
also determine the types of heavy excitations involved in generating of which there will 
be in general several possibilities. Below we will also provide the tree-level diagrams and the 
list of heavy excitations that can generate the leading LNV operators involving two leptons 
(for a systematic study of the dimension 5 operator at one loop see [128]). As in the text, 
in listing the heavy particles we denote by X] , Wj and a heavy vector, fermion 
or scalar with isospin / and hypercharge F, respectively. When the heavy particles can be 
either a heavy vector or heavy scalar with the same isospin and hypercharge, we use B] 
to denote both possibilities. 



A. LNV operators with no quarks 

All the operators below can be generated at tree level, and it is a simple exercise to 
determine the types of heavy physics that can do so. The number of possibilities, however, 
increases rapidly with the dimension of the operator so we will restrict ourselves to those 
with the lowest dimension operators within each group. Still it is useful to note the following: 
if O differs from O' by the presence of 2 derivatives, O ~ D^O', then if O is generated at 
tree level in a certain model, O' will be also generated at tree level. Thus, for example, if a 
model generates a low-energy effective vertex ~ e^W"^ (from D"^ i"^ cfp') , it will also generate 
a Majorana-mass term v\ at tree level. In general, however, it will not generate a e'j^W'^ 
vertex at tree level. As mentioned repeatedly, this particular example is of interest when 
studying Oz//3/3 decay. 

Below we list the LNV operators of dimension < 9 not involving quarks together with the 
sets of heavy excitations that may generate them at tree level. The operators involving two 
leptons can be grouped in 3 sets according to the chirality of the light leptons (operators 
with > 2 leptons will be listed at the end, and are provided for completeness only): 
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a. Two LH leptons plus bosons 



Dim. 


Operator(s) 


5 




7 


£202^ 


9 


£202(]2 



(45) 



Where Q denotes either or 



n 



(^2,00). (46) 

So that, for example, £202^ corresponds to i'^(j)^(f) or D'^l'^cfP'. As mentioned above, each 
entry represents a series of operators obtained by making all possible index contractions 
and having the derivatives operate on all fields. Considering all possible contractions the 
operators of dimension 5 are simply 

The operators of dimension 7 in this category are 



(47) 



^ijab 

;(7-II) 

ijab 

^(7-111) 



^la ^ ^ ]b 1 

Jab ■ yVij , 



^'JLLU -J ' 

= NlCd,^^, , (48) 

as well as C^jjabl^l^O- We will not consider operators of dimension 9 in this category. 

The heavy excitations that can generate (47) and (48) at tree level can be read off Table 
III and are listed in (33). One may think that (48) are still relevant in that they involve 
also the gauge bosons, but in actual calculations any amplitude involving these effective 
operators will have a counterpart involving (47) together with SM vertices. The second 
amplitude will be suppressed by only 1/A, compared to for the first one, with the 

same scale in both cases. As a result all effects of (48) are subdominant in any process; 
accordingly, we ignore these operators in the following. The same holds for operators of the 
form O^^^cp'^ (p) , and for operators of dimension 9 in this category. 

b. One RH and one LH lepton plus bosons 



Dim. 


Operator(s) 


7 

9 


De£03 



(49) 
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a 


c 






/ \ 

b 


d 


1 


a 


b 


c 


d 




e 


e 








e 




e 





TABLE III: Diagrams generating O^^^ at tree level; fields in smaller font inside brackets correspond 
to models with more than one light scalar isodoublet. 

The leading operator of this class has dimension 7 and is given by 

<^Uab = (e^l^Nkb) = ix,XjXkmWe-a W~PRi^b + -- - , (50) 
where the ellipsis denote terms involving the physical scalars and 

=vxi, = 1 , = ^g^v^ (51) 

i 

with g the SU{2) gauge-coupling constant. 

In order to determine the types of heavy excitations that generate O^"^^ we proceed as 
follows. We first write (omitting subindices for simplicity) 

= (^eY<P^t) (<j)^d^4>) +■■■ , (52) 

where the ellipsis denote terms with vector bosons. This term will be generated by graphs 
containing two charged leptons of opposite chiralities and 3 light scalars as external legs. 
A model that generates this graph at tree level must also generate the full operator (at 
tree level) due to gauge invariance. In practice this is obtained by attaching an appropriate 
number of light gauge boson lines to the internal heavy propagators at all points allowed by 
the quantum numbers. The graphs that can generate the term (52) at tree level can have 
two topologies, but chirality prevents one of them; the remaining tree-level graphs are given 
in Table IV and listed in (34). 



33 



a c d 
b e 


1 


2 


a 


b 


c 


d 


e 


]^(3/2) 
^1/2 




£ 


e 




</> 


</> 


x-(3/2) 
^1/2 




e 


e 




<P 


4> 






e 




e 


<P 


4> 






£ 


^ 


e 


<P 


4> 






i 




e 


<P 


4> 






e 


4> 


I 


<P 


4> 


^(1/2) 
^1/2 




e 




e 




4> 


^(1/2) 
^1/2 




e 








i 



TABLE IV: Diagrams generating O^"^^ at tree level; fields in smaller font inside brackets correspond 
to models with more than one light scalar isodoublet. 

c. Two RH leptons plus bosons 



Dim. 


Operator(s) 


7 

9 





The operators e^0^ and e^cj/'cf) vanish for the case of a single scalar doublet. When there are 
more scalar isodoublets these operators generate vertices of the form eeH^H^ multiplied 
by neutral scalars and/or vacuum expectation values. Barring the presence of light single- 
charged scalars the leading operator of this class then has dimension 9 and is given by 

0%, = eaetW^^ ■ Wm • (54) 

The tree-level graphs that can generate this operator are obtained in the same way as above. 
We expand 

0(9) = (ee^) (^(P^dy) {<P^d^4>) +■■■ , (55) 
where the ellipsis denote terms with vector bosons and then look for graphs with two RH 
electrons and 4 scalars in the external lines, which generate this term in O^^); the rest of 
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</> 


</> 
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</> 




T ('1/2) 
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</> 
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T fO) r C(nl 


(1) r m1 
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</> 
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TABLE V: Diagrams generating O^^^ at tree level with three virtual particles; fields in smaller font 
inside brackets correspond to models with more than one light scalar isodoublet. 

the operator will necessarily be generated because we assume the underlying theory respects 
the SM gauge symmetry. There are three diagram topologies presented in Tables V, VI and 
VII. The sets of heavy excitations that can generate the dimension 9 operators at tree level 
are hsted in (35). 

d. More than 2 leptons and bosons 



Dim. 


Operator(s) 


7 








9 









where the middle term in the last line represents objects of the form Di'^i(l)'^, i'^ie(l)'^, etc. 

The important conclusion of the above arguments is that there are only three interesting 
LNV operators of lowest dimension: 0^^\ 0^'^\ 0'^'^\ Aside form their LNV effects these 
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TABLE VI: As table V but with a different topology. 
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TABLE VIL As table V but with only two virtual particles, 
operators in general transform non-trivially under CP. 
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